Quantum phases of light coupled to a periodic Bose-Einstein condensate in a cavity 
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We study a composite photon-Bose condensed atoms (BEC) system in an optical lattice confined 
in a cavity. We show that the system is able to produce an effective photon-photon repulsion in the 
presence of an external pump. We predict for the first time zero-temperature quantum phases of 
light coupled to BEC in the photon blockade regime in the presence of an external pump. 
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The search for interesting and potentially useful 
quantum-mechanical phenomena on a mcsoscopic scale 
in condensed matter and atomic physics is a challenging 
task. The combination of cold atoms and cavity QED 
(quantum electrodynamics) is a conceptually new regime 
of cavity QED, in which all atoms occupy a single mode 
of a matter-wave field and couple identically to the cav- 
ity light field, sharing a single excitation. Experimen- 
tal implementation of such systems has made significant 
progress [l|, 0, The atomic back action on the field 
introduces atom-field entanglement which modifies the 
associated quantum phase transition 0]. Scattering of 
light from different atomic quantum states creates dif- 
ferent quantum states of the scattered light, which can 
be distinguished by measuring the photon statistics of 
the transmitted light. [!, @, 0, H- The band struc- 
ture of the intracavity light field has been shown to in- 
fluence the structural properties of the BEC [9(. Ex- 
otic phases of ultracold atoms in cavities have been pre- 
dicted recently [10] . The optical nonlinearities in the 
Jaynes-Cummings model generated due to the coupling 
between the atom and the photons leads to an effective 
photon-photon repulsion [l! IE!. However it was shown 
that photon-photon repulsion degrades in the presence of 
many atoms 14|. Adding photons to a two-dimensional 
array of coupled optical cavities each containing a single 
two-level atom in the photon-blockade regime, a long- 
lived, strongly interacting dressed states of excitations 
(coupled atom-photons) are formed which can undergo 
at zero temperature a characteristics Mott insulator to 
superfluid quantum phase transition 15]. Here we de- 
scribe a coupled photon-atom system formed by adding 
photons to an one-dimensional array of Bose condensed 
atoms (BEC) confined in a cavity in the photon-blockade 
regime that exhibits photon-photon repulsion in the pres- 
ence of an external pump. We search for a Hubbard- 
model type interactions in the coupled photon-atom sys- 
tem and predict the existence of quantum phases of com- 
posite photon-atom states at zero temperature. 
We consider an elongated cigar shaped Bose-Einstein 
condensate of N two-level 87 Rb atoms in the \F = 1 > 
state with mass m and frequency uo a of the \F = \ >—t 
\F' = 2 > transition of the Z?2 line of s7 Rb, strongly 




BEC in 
optical lattice 



Transmitted 
light 



FIG. 1: Setup of the system. Atoms in a lattice are illu- 
minated by the pump and the scattered light is transmitted 
through one of the mirrors and measured by a detector. 



interacting with a quantized single standing wave cavity 
mode of frequency lu c . The internal cavity field is linked 
with the input and the output can then be calculated 
using the boundary condition at the cavity mirror. In 
order to create an elongated BEC, the frequency of the 
harmonic trap along the transverse direction should be 
much larger than one in the axial (along the direction of 
the optical lattice) direction. The system is also coher- 
ently driven by a laser field with frequency uj p through 
the cavity mirror with amplitude r\. The two sided cav- 
ity has two partially transparent mirrors with associated 
loss coefficient 71 and 72 . The cavity decay is assumed to 
dominate over the spontaneous decay of all atoms thus al- 
lowing us to omit the effect of atomic decay. This system 
is modeled by a Bose- Hubbard (BH) type Hamiltonian in 
a rotating wave and dipole approximation Adiabati- 
cally eliminating the excited state and retaining only the 
lowest band with nearest neighbor interaction, we obtain: 
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J d 3 xw(f — fj) cos 2 (kx)w(r — ^j±i) and w(r — rj) is the 
Wannier function. The BH Hamiltonian of Eqn. (1) is 
different from the one derived in ref. [4|, where the cav- 
ity photon operator is not dependent on the lattice site. 
A recent experiment demonstrating that BEC atoms 
could be loaded into a single site of an intracavity optical 
lattice and deterministic transfer of the BEC atoms into 
successively addressed single sites of the lattice, each of 
which is differently coupled to the cavity photons pro- 
vides a justification for the site dependence of the cavity 
photon operator. A single photon can be pumped into a 
well of the optical lattice by a low intensity pump laser 
with width less than the lattice spacing so that the spa- 
tial overlap of the pump and the lattice is no larger than 
the lattice spacing. The first photon will form a compos- 
ite photon-atom system and a secong photon that comes 
in will experience a repulsion unless there is a second 
atom to form a second photon-atom system. In ref. [l6| , 
a probe laser with 3.9/im width and lattice spacing of 
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6Afim was utilized. The parameter Uq — -— is the op- 

tical lattice barrier height per photon and represents the 
atomic backaction on the field Here we will always 
take Uq > 0. In this case the condensate is attracted to 
the nodes of the light field and hence the lowest bound 
state is localized at these positions. Along x, the cavity 
field forms an optical lattice potential of period A/2 and 
depth hUo(ajCij) at each lattice site. Here A a — lo v — to a 
and A c = lu p — uj c are the large atom-pump and cavity- 
pump detuning, respectively and A c > A a . In this work 
we will consider only the case A a > 0. The atom-field 
coupling is written as g(x) — gocos(kx). Here dj and 
bj are the annihilation operator for a cavity photon and 
the bosonic atom in the j th site of the optical lattice re- 
spectively. Here a s is the two body s-wave scattering 
length. The onsitc energies Jo and Eq are set to zero. 
The Heisenberg equation of motion for the bosonic field 
operator b is 



bj = -iU (l + a)a^ J \b J+ i + j - y |S J+ i + bj-ij 

The behaviour of the internal cavity mode is obtained 
from the quantum-Langevin equation which for a single- 
mode cavity becomes 



set of coupled equations describing the dynamics of the 
compound system formed by the condensate and the op- 
tical cavity. We will work in the bad cavity limit, where 
typically, *yi and 72 are the fastest time scale ( this means 
that the cavity decay rates are much larger than the os- 
cillation frequency of bound atoms in the optical lattice 
of the cavity ). In this limit the intracavity field adia- 
batically follows the condensate wavefunction, and hence 
we can put dj — 0. We treat the BEC within the mean 
field framework and assume the tight binding approxima- 
tion where we replace bj by <pj and look for solutions in 
the form of Bloch waves (f>j = Ukexp(ikjd)exp(—ifit/h). 
Here /i is the chemical potential, d is the periodicity of the 
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lattice and — ^7 bjbj = \uk\ = no (atomic number den- 
sity) , / is the total number of lattice sites and J2j n o = N 
(total number of atoms). The relationship between the 
input and output modes may be found from using the 
boundary conditions at each mirror, d ou t(to) + di n (uj) = 
yfyiaj(u) and b out {uj) + b in (u>) = ^72%' (w) [ll|. In fre- 
quency space, we obtain 



7i?? + (7i/2 - 72/2 + iA)a in (uj) + y/ji 72 b in (w) 



Where d out ( in) (uj) = 
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a out ( in ){t)dt and 



A = A c + u> — 2 JuqIIq cos(kd) We find that due to 
the atomic backaction, the quantum state of the out- 
put cavity field develops a band structure due to the 
strong coupling with the condensate, analogous to pho- 
tonic band gap materials 17| . The average photon num- 
ber < al ut d ut > measures the light transmission spec- 
tra and is different for the Mott insulator (MI) and the 
superfluid phase (SF) Q. If the two mirrors are the 
same, 71 = 72 = 7 and near resonance A « 0, we 



find: d out (uj) 



The resonance point 



{ 7 -i A'} 

is one where the superfluid fraction is minimum [9( ■ This 
shows that the cavity now behaves like a shifted through- 
pass Lorentzian filter. The input field will be completely 
reflected if A >> 7 (when the atoms are in the deep 
superfluid regime), a ou t(uj) ~ — hi n (uj). The state of the 
BEC together with the cavity parameters control the out- 
put field of the cavity. 

The Hamiltonian of Eqn.(l) in the mean field of the 
atoms is written as 
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H MF = 2En cos kd + 2hU Q Jn Q cos kd(l + a\aj) 
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Here &i n and bi n are the external input fields incident 
from the two mirrors.. Equation (3) and (4) represents a 



Analogous to the Jaynes-Cummings model, the mean 
field Hamiltonian H MF leads to photon-photon repulsion 
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FIG. 2: Dimensionless eigenspectrum, (E± — E c )/hUo as a 
function of A/Uo for n = 1. The eigenspectrum splits into two 
branches , corresponding to the dressed states, E+ t „ (upper 
branch) and E- ]Tl (lower branch). The two branches anti-cross 
at A/Uo = 0, with the splitting increasing with increasing 
pump strength. 
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FIG. 3: Dimensionless eigenspectrum, (E± — E c — (2n — 
l)A)/hUo versus A/Uo for a fixed pump strength and n = 
1 (solid lines) and n = 2(dashed lines). The two branches 
anti-cross at A/Uo = 0, with the splitting increasing with 
increasing excitation number, n. 



in the presence of the pump. The conserved particles in 
out model are not the photons but composite photon- 
atom states. The eigenvalues of the H^ F (mean-field 
Hamiltonian at each lattice site)in the basis \g, n > and 
\g, n — 1 > is found as, 



E±, n = {2n 2 1)HA ± \Vh 2 ^ 2 + 4rcfiV + f (6) 

where \g,n > is the ground state of the compos- 
ite atom-photon system with n excitations. Here 
A = 2UquqJ cos kd — A c and E c = 2uq(E cos kd + 
hUo Jcos kd) + U/2\uk\ 2 . In Fig. 2, we show the dimen- 
sionless eigenenergies (E± — E c )/HUq as a function of 
A/C/o for n = 1. The bold lines are for n/Uo = 1 and 
the dashed lines are for rj/Uo = 2. The increasing en- 
ergy separation with increasing pump strength is an in- 
dication of photon-photon repulsion. In Fig. 3, a plot of 
(E± — E c — (2n— l)A)/HUo versus A/C/o for a fixed pump 
strength and n — 1 (solid lines) and n = 2 (dashed lines). 



The on-site photonic repulsion is evinced by the increas- 
ing energy separation with n. In the absence of the pump, 
the photon-photon replusion is absent. 

The effective Hamiltonian for our extended Hubbard- 
like system is given by a combination of the mean-field 
Hamiltonian with photon hopping between wells and the 
chemical potential term. We introduce a superfluid order 
parameter tp —< dj >, which we take to be real and use 
the decoupling approximation, aldj =< fit > dj+ < 

a,j > cij— < fit >< cij >. The resulting effective mean- 
field Hamiltonian can be written as a sum over single 
sites 



ttMF 
H eff 
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where t is the nearest neighbor hopping energy of the 
photons coupled with the atoms (actually it is the exci- 
tations which hop from one well to the other) and jij is 
the chemical potential at site j. We assume zero disor- 
der such that fXj = fi for all sites. Here z = 2 is the 
number of nearest neighbors. To obtain the system's 
zero-temperature phase diagram, we use the procedure 
of ref. [l8(. This Hamiltonian is diagonal with respect 
to the site index j, so we can use an effective onsite 
Hamiltonian. If we introduce Uq = Uo/zt, fi = fb/zt, 
Hf IF = Hf IF /zt, we find H5f = H° + ipV with H° = 

j + cij). The un- 
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, = Hf IF /zt, we find H*}f = H° 
MF ' ^ 2 fi{a]a 3 + n Q ) and V = -(fi} 



perturbed ground state energy of the state with exactly 
n particles is E® n = E_ <n — fin, where E_ <n = E_ <n /zt. 
We only need to consider the negative branch for the pur- 
pose of determining the ground state since < E + , H . 
A change in the total number of excitations per site will 
occur when E— tn +i — Ji{n + 1) = E_ >n — fin. We can de- 
termine the critical chemical potential, /2 c (n), where the 
system will change from n to n + 1 excitations per site as 
fi c (n) = A+ (xn-Xn+i), where Xn = v/A 2 /4 + nh 2 rf / zt. 

Since the interaction V couples only to states with one 
more or less excitations than in the ground state, we find 
the second order correction to energy 



#(2) = y_ 

9 A + (x g - Xg-x) - A* 
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-A + (Xg+l -Xg) + V 
(8) 

According to the Landau procedure for second-order 
phase transitions, we write the ground state as an expan- 
sion in V as E g {ip) = a (g, p.,fj, A) + a 2 (g, fl,fj, A)ip 2 + 
0(ip^). E g (ip) is minimized as a function of the super- 
fluid order parameter ip. We find that ip = when 
02(3, fi, fj, A) > and that if) ^ when ai{g, fi, fj, A) < 0. 
This means that a 2 {g, fi, fj, A) = signifies the boundary 
between the superfluid and insulator phases of light. This 
yields 
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FIG. 4: Phase diagram of the coupled periodic photon-BEC 
system in the cavity in the presence of the pump. The ver- 
tical axis shows fi/zt and the horizontal axis shows f\. The 
plot is shown for g = 1,2 and A = —1,0,1. Dominating 
the right-hand edge (where photonic repulsion dominates at 
higher pump strength) is the Mott-insulator (denoted by MI) 
and the superfiuid phase is found on the left-hand edge (de- 
noted by SF). 



ber should be below 7 2 /2<?o ~ 0-04 [HI- Photon loss can 
be minimized by using high-Q cavities and thus ensur- 
ing that the light field remains quantum-mechanical for 
the duration of the experiment. It is important that the 
characteristic time-scales of coherent dynamics are signif- 
icantly faster than those associated with losses (the decay 
rate (7) of state-of-art optical cavities is typically 17 kHz 
[2fj|). In summary, we have shown that a system of peri- 
odic Bose condensed atoms coupled to cavity photons in 
the photon blockade regime exhibits a rich strongly corre- 
lated dynamics. The composite photon-atom states thus 
formed undergoes a characteristic Mott insulator to su- 
perfiuid quantum phase transition. Because of the Mott 
phase's robustness, the present system provides a useful 
platform for developing concepts in quantum information 
processing. 



/2± = i(2A-x 9 +i+X 9 -i-l)±^l + + 2(2n+ 1)E, 

(9) 

where £ = Xg+i+Xg-i~ %Xg- The subscript ± denotes 
the upper and lower halves of the superfiuid regions of 
phase space. Fig. 3 shows a plot of Eqn.(27) for g = 1, 2 
and A = —1,0,1 as a function of fj. By equating /2_ 
and /2_|_ we can find the point of largest fj max for each su- 
perfiuid region (SF). The dynamics illustrated in Fig. 3 is 
extremely rich. Dominating the right-hand edge (where 
photonic repulsion dominates at higher pump strength) 
is the Mott-insulator (denoted by MI) and the super- 
fluid phase is found on the left-hand edge (denoted by 
SF).The size of the superfiuid region is found to increase 
with the excitations. On increasing the pump strength, 
the optical potential increases and hence the atoms loose 
their superfiuid behavior Q. The transition from the 
SF to the MI phase of the excitations occurs at 7] ma x- 
In ref. 0] it was shown that the superfiuid fraction of 
the cold atoms shows a minimum at A = 0. In accor- 
dance with this result, we find that the size of the su- 
perfiuid lobes of the composite photon-atom system is 
minimum at A = 0. To find the eigenenergies and to 
experimentally identify the various phases of the coupled 
BEC-cavity system, one can perform a transmission spec- 
troscopy with the scattered light by direct read out of 
the number of photons coming out of the cavity through 
Eqn.(6). The transmission of the scattered light is mon- 
itored as a function of A. To probe the system in the 
weak excitation limit, the mean intracavity photon num- 
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